Geometric phase for an adiabatically evolving open quantum system 
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We derive an elegant solution for a two-level system evolving adiabatically under the influence of 
a driving field with a time-dependent phase, which includes open system effects such as dephasing 
and spontaneous emission. This solution, which is obtained by working in the representation cor- 
responding to the eigenstates of the time-dependent Hermitian Hamiltonian, enables the dynamic 
and geometric phases of the evolving density matrix to be separated and relatively easily calculated. 

PACS numbers: 



I. INTRODUCTION 

The discovery by Berry 0, Q that a (non-degenerate) 
state of a quantum system can acquire a phase of purely 
geometric origin when the Hamiltonian of the system im- 
dergoes a cyclic, adiabatic change has lead to an explo- 
sion of interest in this and related geometric phases in 
quantum mechanics, both from a theoretical perspective, 
and from the point of view of possible applications, the 
latter including applications to optics (where the geo- 
metric phase was first discovered jj]), NMR and molec- 
ular physics, and to quantum computing |j, Q- Since 
Berry's work, and the demonstration that Berry's phase 
can be understood as a holonomy associated with the 
parallel transport of the quantum state [3, there have 
been numerous proposals for generalizations. The first 
of these was due to Wilczek and Zee p who, by consid- 
ering a Hamiltonian with non-degenerate eigenstates, es- 
tablished the existence of an intimate connection between 
Berry's phase and non-Abelian gauge theories. The re- 
striction to changes occurring adiabatically was relaxed 
in the work of Aharonov and Anandan JT"! while Anandan 
Q generalized the geometric phase to the non-adiabatic 
non-Abelian case. The restriction of cyclicity was re- 
moved by Samuel and Bhandari ^ and by Pati [Ij . All 
of this work is concerned with geometric phases of pure 
states of closed systems and is now standard, though JJ| 
indicated extensions to taking account of quantum mea- 
surements and consequent non-unitary evolution. A nice 
overview, theoretical as well as experimental, is given in 

M- 

More recently attention has turned to studying geo- 
metric phases for mixed states, though there is not yet 
a standard description for geometric phases associated 
with mixed states. 

As realistic systems always interact with their environ- 
ment, and as an open system is almost always to be found 
in a mixed state, open systems are a natural source of 
problems involving the geometric phases of mixed states. 
Garrison and Wright [ij were the first to touch on this 
issue in a phenomenological way, by describing open sys- 
tem evolution in terms of a non-Hermitean Hamiltonian. 



This was, in fact, a pure state analysis, so it did not, 
strictly speaking, directly address the problem of geomet- 
ric phases for a mixed state, but this work raised issues 
which could potentially have a bearing on the analysis of 
the mixed state problem. In fact, they did point out that 
a proper treatment of an open system would require mak- 
ing use of the density operator approach. Nevertheless, 
they arrived at an interesting result, a complex geomet- 
ric phase for dissipative evolution. This is a result that 
has been recently put into doubt by a master equation 
treatment by Fonseca Romero et al and Aguiar Pinto and 
Thomaz 13, 14]. 

The first complete open systems analyses of geometric 
phase for a mixed state, from two different perspectives, 
is to be found in the papers of EUinas et al jl5j| and Gam- 
liel and Freed [l6| . The former worked with the standard 
master equation for the density operator of a multilevel 
atom subject to radiative damping and driven by a laser 
field with a time-dependent phase. What is of interest 
in their approach is that it entailed introducing eigenma- 
trices of the Liouville superoperator of the master equa- 
tion for the damped system. The system Hamiltonian 
was allowed to vary adiabatically, with the result that a 
non-degenerate eigenmatrix acquires a geometric phase 
as well as a dynamic phase. In [lal , the effect of the envi- 
ronment was modelled as an external classical stochastic 
influence which, when averaged, gives rise to the relax- 
ation terms of the master equation for the system. In 
both cases the effects of any geometric phase was then 
shown to be present in measurable quantities such as the 
inversion of a two state system. 

Since then, research has been increasing rapidly into 
the problem of defining a geometric phase for mixed 
quantum states for both unitary and non-unitary evo- 
lution, motivated to a very large extent by the need to 
understand the effects of decoherence in quantum com- 
putational processes that exploit geometric phases as a 
means of constructing intrinsically fault-tolerant quan- 
tum logic gates. This issue has been addressed from two 
points of view, the first holistic in nature wherein the 
aim is to identify a geometric phase to be associated with 
the mixed state itself, and the second, essentially the ap- 



proach of |l3 , which works with the pure state geometric 
phases of an appropriate set of parallel transported basis 
states, which then gives rise to geometric phase factors 
in the off-diagonal elements of the density operator. No 
geometric phase is explicitly associated with the mixed 
state itself; instead observable quantities that will exhibit 
the effects of the geometric phases of the underlying basis 
states are determined. 

The former, holistic approach was first introduced in a 
formal way by Uhlmann !l7|, and in a different way, based 
on phase-sensitive measurements via interferometry, by 
Sjoqvist et al |l3| for unitary evolution of a mixed state, 
and later for non- unitary evolution |19l |20| . The phase 
defined in this way is not the same in all respects to that 
proposed by Uhlmann ^^ . 

The latter kind of approach has been used only for 
open systems, and involves working with the (Markovian) 
master equation of the open system. The approaches 
used involve either solving the master equation of the sys- 
tem [l3l llil IIEI llH ^ or employing a quantum trajectory 
analysis |22l l23l |24 || to unravel the dynamics into pure 
state trajectories, and calculating the geometric phases 
associated with individual pure state trajectories. Noise 
of a classical origin, such as stochastic fluctuations of the 
parameters of the Hamiltonian have also been studied by 
Chiara and Palma [23 • In essence, the common feature 
is not so much to propose a new definition of geometric 
phase for a mixed state as to show how the underlying 
existence of a geometric phase will nevertheless show up 
in the observed behaviour of an open quantum system. 
It is this perspective that is adopted in the work to be 
presented in this paper. 

Here we introduce an elegant approach for solving the 
central master equation which is based on introducing a 
unitary transformation due to Kato pq described in the 
classic text by Messiah j27j . 

The new picture is defined via a time-dependent 
unitary transformation v4^(i), usually referred to as a 
rotating axis transformation, which is such that the 
transformed system Hamiltonian has time independent 
eigenspaces. This method is extended by showing that 
under the conditions of adiabatic evolution, all the infor- 
mation on geometric phase for a closed loop is contained 
within A{t), and is regained by transforming back to the 
original picture. The goals of this approach is its sim- 
plicity, since one needs only to calculate the geometric 
phases for the eigenstates of the Hamiltonian H{t), and 
the fact that dynamic and geometric phase are separated 
in a clear way. In fact, perfect separation of geomet- 
ric and dynamical contributions is obtained provided the 
Hamiltonian evolution is adiabatic and the coupling to 
the environment is weak. This approach bears some sim- 
ilarity to that used by Fonseca RomeroFonseca^lSj who 
make use of several unitary transformations to separate 
the geometric phase from the dynamic phase. However, 
the rationale for their transformations, and the origin of 
the geometric phase, is somewhat elusive in their anal- 
ysis. In contrast to |l3| . with the transformation intro- 



duced here, the parallel transport condition is essential 
and explains the appearance of the geometric phase. 

Within the approach used here, it is possible to show 
explictly how to achieve, under certain circumstances, 
operational cancellation of the dynamic phase, thereby 
making the geometric phase accessible in experiments. 
An example of where this is possible is given in Sec- 
tion |IlTl 

This paper is organized as follows. In Section ^ we 
present the main ideas. In Section IIIII we look at sev- 
eral examples. In Section we summarize our results 
while in Section llVI we indicate possible new directions, 
including generalisations to non-cyclic evolution and non- 
Abelian holonomies. An analog to the adiabatic theorem 
is proved for a general Lindblad equation 28\ in the Ap- 
pendix. 



II. THE ROTATING AXIS TRANSFORMATION 
FOR NONDEGENERATE MULTILEVEL 

SYSTEMS 

A. Geometric Phases for a Closed System 

For the present we consider the case of a closed system 
so as to introduce the basic method employed here. Sup- 
pose we have a system with Hamiltonian H{t), a function 
of time due to the dependence of H on parameters whose 
values can be changed in time. This Hamiltonian will 
have instantaneous eigenvectors \n(t)) with eigenvalues 
Enit): 



H{t)\n{t)) ^En{t)\n{t)) 



(1) 



For simplicity we assume Ei{t) 7^ Ej{t) for i ^ j. This 
restriction will be removed in section IV to obtain non- 
Abelian holonomies. For an adiabatically slow change 
in the system parameters, these eigenvectors will also 
change in such as way as to satisfy the parallel trans- 
port condition Q]: 



{n{t)\j^\n{t)) 



0. 



(2) 



At this point we introduce a unitary operator A{t) via 
the equation 



AiDHO)) ^\n{t)) 



(3) 



This completely defines A{t). Note that because of the 
path dependence of the parallel transported eigenstates 
\n{t)) , the operator A{t) has a non-integrable nature, 
and, as we see later, will contain the information on the 
geometric phase. 

This unitary operator can now be used to remove the 
time dependence of the eigenstates of the Hamiltonian. 
Thus, if we define 



H"^ = A^HA 



(4) 



we note that the eigenvectors of H^ 
and hence are time independent. 
Schrodinger equation is then 



are now just |n(0)) 
The transformed 



for 



Ct[p\ = 2rpri' - r^Tp - pr+r 



(11) 



H^\^^)^^n(A^A\^^) + ^\i,^) 



dV 



(5) 



where |'0'^) = ^4^1?/;). If \n{t)) is parallel transported, 
then, to the lowest-order adiabatic approximation, one 
can neglect the terms containing A\t)A{t) so that the 
transformed Schrodinger equation becomes |27| 



H^W 



m^y) 



(6) 



The solution of Eq. ^ contains no geometric contribu- 
tion — it gives the dynamic contribution to the phase of 
any adiabatically evolving state. So we have extracted 
the dynamics from the geometric influence of the time- 
varying Hamiltonian. The geometric contribution is en- 
tirely contained within the operator A. To obtain this 
information we have to transform back to the original 
picture in terms of the states \'ip): 



IV') = A\i;' 



(7) 



If the Hamiltonian undergoes a closed loop in time T 
i.e. H{T) = H{0), then the parallel transported eigen- 
states \n{T)) return to the initial eigenstates |n(0)) 
up to the geometric phase. Hence we have A{T) 
diag(e''^i ,•■•,< 



and where p{t) is the density operator for the system of 
interest, p{t) is its derivative with respect to time, H(t) 
is the system Hamiltonian and the dissipation operators 
Fq arise due to the presence of the reservoir. To obtain a 
geometric phase we once again consider an adiabatically 
changing Hamiltonian H{t), with H{T) = H{Q), as in the 
preceding section. Upon introducing the operator A{t) 
defined in Eq. (jJJ we can transform the master equation 
with the unitary operator At(t), which leads us to the 
new master equation 

/ (t) =^[H^{t), p^{t)] + p^{t)A\t)A{t) 
in 

k 

^At(i)i(i)/(i) + i^£p.[p^(f)]. (12) 

q:=1 

Note that F^ = A^TaA is in general time-dependent even 
when Fq, is not. As shown in the Appendix, the terms 
containing A\t)A{t) can be neglected, as in the case of 
a unitary evolution but with the additional requirement 
that the damping is weak. This establishes the parallel 
transport condition Eq. Q for weakly damped systems 
in an adiabatic evolution. We then get 



) where ^„ is the geometric phase as- ^A(^) ^ L[H^{t),p^{t)] + i V £r^[/9^(i)]. (13) 



sociated with the eigenstate \n{T)) . This result can be 
readily generalized if the system is in a mixed state p 
Introducing the notation 



p^{t) = A\t)p{t)A{t) 



(8) 



we obtain 



p{T) ^A{T)p^{T)A\T) 



iipi 



=diag(e 
X diag(e" 



„1<PJVN -4 



)p^(r) 

')■ 



(9) 



This holonomic transformation multiplies the off- 
diagonal elements of the density operator p^ by a phase 
g^iivj-Vi) ^ which is the difference of the geometric phases 
of the eigenstates of the Hamiltonian H{T). 



ifi 

Q = i 

The solution of Eq. ^VA\ contains no geometric contri- 
butions. To regain p{T) we have to transform back to 
the original picture as in the unitary case, Eq. JHJ. 

p{T) =A{T)p^{T)A\T) 

=diag(e''^S...,e''^")/(r) 

xdiag(e-''^S...,e~*'^") (14) 

with Lpn being the geometric phase for the state \n) . 

III. EXAMPLES FOR TWO-LEVEL SYSTEMS 
A. Optical Resonance with Spontaneous Emission 



B. Open System and Master Equation 

Systems that are coupled to a reservoir (or environ- 
ment) can usually be described by a reduced density 
operator that evolves according to a master equation, 
which, in many cases, can be written in the Lindblad 
form ,28,] 



m 



ih 



[H{t),p{t)] 



1 '^ 



,[p(t)l (10) 



We consider a two level atom in a classical reso- 
nant laser field. In the rotating-wave approximation the 
Hamiltonian for this system is 



H 






(15) 



The detuning A, the coupling strength 57 and the phase 
4> are properties of the laser. To induce a geometric phase 
we change the phase 4>{t) of the laser field slowly in com- 
parison to E/h = (fi^ -I- iA^)2j which is the absolute 



value of the eigenenergies of the Hamiltonian divided by 
h. The eigenvalue equation is 



H{t)\ + {t)) =E\ + {t)) 
Hit)\-{t)) ^-E\-{t)) 



(16) 
(17) 



with 



\ + {t)) ^e-^^f'^'^'^'^cos-le) +e 



»0(*)cos^ I singly) 



|-(t)) ^_e-»-^(*)™-'lsin-|e) + e^-^^^'"' i cos -jg) 



(18) 

\9) 
(19) 



and 



sm - = 
2 



cos ■ 



'E-^hA 



2E 



' E- 



\hA 



2E 



(20) 



(21) 



|e) and \g) denotes the excited state and the ground state 
of the two level atom, respectively. Note that |+) and 
|— ) satisfy the parallel transport condition as required in 
Section ini 

Furthermore wc want to include spontaneous emission 
as a source of dissipation. In the weak coupling limit the 
master equation is known to be 



p{t)^^[H{t),p{t)] + \cT[p{t)]. 



for 



\f\(T- =\/A 





1 



(22) 



(23) 



Here A denotes the spontaneous emission rate. The task 
here is to solve Eq. H22|l in the adiabatic and weak damp- 
ing limit. As in Section |n] we define the operator A{t) 
by A{t)\ ± (0)) = I ± (i)) . After the transformation of 
Eq. H22|l with A^(i), the Hamiltonian is not diagonal. 
Hence we carry out another transformation with an op- 
erator B'^ that is defined by 



i?t| + (0))=|e) 



(24) 



As i?^ is time- independent, we obtain no term B'^ B in the 
master equation. We can carry out both transformations 
together with the operator C^(i) = B'^ A^it), which turns 
out to be 



for 



and 



H^ = C\t)H{t)C{t) = 



E 
--E 



(27) 



T^{t) =CHt)TC{t) 



-^V~x( 



\ cos^ ^Q^'t'W '^°^ ^ 



sin'^ ^g-#(t) cos( 



(28) 



We show in the Appendix that only the absolute value 
of the off-diagonal elements of the dissipation operator 
r'-^ (t) contributes in the solution of Eq. H26|) . Using only 
the absolute values of these elements does not simplify 
the calculation, but multiplying them with a time inde- 
pendent factor e^^ and averaging over all < /? < 27r 
does. Doing so, we first rewrite Eq. (|26|l 

In Eq. H29() we substitute the Lindblad operator T^ with 
the new Lindblad operators F? with 



-pL. _ ./ \ ( -"" 2 " 2 " 2 



2 9_„-if3 



c _ p-lcos^sm^ sm" ^e 



\ cos^ ^e ' 



< /3 < 27r 

(30) 

as explained previously. Finally we get a new master 
equation, which is equivalent to Eq. H2()|) in the adiabatic 
and weak damping limit: 



pCit) ^l_^H^^p^{t)] + ^j^\rc^,y{t)]dl3 (31) 

Eq. (|31l) may seem to be more complicated than Eq. I|26|l . 
but if we evaluate the term behind the integral in Eq. H31|l 
we realize that many contributions are cancelled out by 
the integration. If we write 



then Eq. (|31|l reduces to the two independent equations 



a{t) = A 



sm a \ sm — h cos - 

2 V 2 2 



m 



2i 



Eb- Xb { sin^ - cos^ 



e 1 

2 + 2 



(33) 

(34) 



with the solutions 



C{t) 



'i(f>{t) sin 



2 COS : 



(t)cos^lgjj^J 



,—i(p{t) cc 
i(p{t) sin 



2 sm ■ 
8 e ■ 

2 COS o 



(25) 



The master equation for p''^' {t) — C^{t)p{t)C{t) is, from 
Eq. (O and Eq. ^, 



P'^it) 



\H^,p^it)] + l-Crc^,)[p^it)] (26) 



ih 



a{t) = a(0) 



sin'^ I -I- cos4 I 



,4 e 



sin'^ I + cos* I 



b{t) =6(0)e-2.iig-At( 



^At(sin''f+cos''f) 



(35) 
(36) 



As the last step we need to evaluate p{t) = 
C{t)p^{t)C\t). As the inversion provides an operational 
quantity for inferring the geometric phase by measuring 
the relative proportion of ground vs excited states, and 
because the terms become rather long, we only write the 
inversion 'w{t), which is 

w{t) =pii - P22 = (2a(i) - l)cos6' 

-2sin6lRc(5(i)e*'^(*)™^'') (37) 

To compare this result with that found by EUinas et al 



|15| . we set /o(0) = 5+^0-3 with \p\ < i and substitute for 
sin I and cos | from Eq. H20|l and Eq. (|21|l . respectively. 
Furthermore we define 



K 



2n^ 



4^2 



and G 



8172 + 2A2' 



(38) 



If we furthermore consider the inversion at a time T at 
the end of the cyclic evolution we finally get for the in- 
version: 



wiT) 



2Kp+l^_j,^T 




2Kp 



f2ET 



\ fi 



2tt e" 

2EJ 



GXT 



~E 



(39) 



which is the same as that derived by |15l |. The dynamic and geometric phases are found in the cosine term in this 
expression: the difference of the dynamic phases of the eigenstates of the Hamiltonian is given by 2ET/h, and the 
difference of the geometric phases of these eigenstates (for 4>{T) = 2tt) given by 27ThA/2E. This term is diminished 
by a damping factor exp(— GAT) which can influence the observability of the geometric phase effect on the inversion. 
The issues of time scales to observe the effect of the geometric phase have been discussed in jl3|. However, for the 
present, we wish to point out that the result above has been derived here by use of a simple transformation into a 
rotating frame. This is to be contrasted with the nruch more complicated approach of [l5J, based on calculating the 
eigenmatrices of the Liouvillian. 



B. Optical Resonance with Dephasing 



As in the previous subsection we treat a two level atom driven by a resonant electromagnetic field. This time we 
assume the damping is due to dephasing that occurs as a consequence of phase changing collisions, which changes the 
relative phase between the excited state and the ground state of the atom (in contrast to strong collisions that change 
populations of eigenstates). Since the phase change can vary for each collision we have to consider a one dimensional 
manifold of dissipation operators 



T^ = ^/X{a)( 



1 

e*" 



— TT < a < TT, 



where A(a) is the dephasing rate density. Hence we get the master equation 

m = l[Hit),p{t)] + r (r„p(<)rl - Xia)pit)) da 



(40) 



(41) 



with the Hamiltonian H{t) from the previous subsection with a slowly changing phase (j)(t) again. Thus, we get the 
same parallel transported eigenstates of the Hamiltonian and we can start by carrying out the same transformation 
as in the previous subsection. In the adiabatic and weak damping limit we obtain the transformed master equation 



p^it) =l.[H^\p^'it)] + r {r^{t)p^it)T^Ht) - X{a)p^it)) da 



(42) 



for 



and 



H^ = C\t)H{t)C{t) = 



T-it) = cHt)T.cit) = vw) ( ,3in7sin .er(,(r):;si.,) 



E 
-E 



cos' 



sm 



isin6'sin|e^(^'^(*)'=°"^+i 
sin^ I + e*" cos2 | 



(43) 



(44) 



Now everything is much the same as in the previous subsection. Finally we find for the components a and b of the 
density operator p'^{t) the decoupled differential equations 



.. 917.9C _. 9(7.9(7 

a = —Afa cos - sm - + 2 f cos - sm - 

2 2 •^ 2 2 

-2E 



gb ( sin - — cos - ) ) — /& ( sin - + cos"* - 

r 



(45) 
(46) 



where 



and 



/ = 



/+7r 
A(a)(l — cos a) do; 



/+TT 
X{a.) sin a do; 
-7r 



(47) 



(48) 



are properties of the model describing the damping col- 
lisions. The solutions of these equations are 



a{t) = (a(0) - i) e-4/tcos^ §-n^ f + 1^ (49) 

b{t) - 6(0)e<"^+s(^"^' i-™^' i))*e-(^"' ^+™^' ^)^* 

(50) 

To calculate the inversion w^t) we can take Eq. H37|) 
and substitute a(t) and 6(i) with Eq. l(^ and Eq. (|5n|l . 
respectively. Using Eq. ipiljl and Eq. (|?T)l as well as the 
previous definition of K we finally find at time t — T: 



w(T) =2p 




X cos 



-KfT 



hA 2E\ hA\ \ , , 



This result is similar to that found in the case of sponta- 
neous emission in Eq. H39() . There appears in the cosine 
term in Eq. (|51|l the difference of the dynamic phases of 
the eigenstates of the Hamiltonian, as given by 2ET/h, 
and the difference of the geometric phases of these eigen- 
states (for 0(T) = 27r), given by 2iThA/2E. This term is 
also diminished by a damping factor exp(—KfT) which 
influences the observability of the geometric phase effect 
on the inversion. Both this term and an additional con- 
tribution of a shift in the Rabi frequency by g^ arise 
through the presence of dephasing, though the latter will 
only appear if the dephasing rate density X{a) is not 
symmetric. 



C. Spin in Magnetic Field with Dephasing 



As another example we consider the simple model of 
a spin-i particle in a magnetic field with constant field 
strength, which demonstrates how to remove the dy- 
namic phase in a standard model. To induce a geomet- 
ric phase we change the direction of the magnetic field 



slowly in comparison to E/h. As a source of decoherence 
we consider dephasing which is defined by 



r„(i)|e(i)) = VWMt)) 



(52) 
(53) 



\g{t)) and \e{t)) are the parallel transported eigenstates 
of the Hamiltonian with spin parallel and antiparallel 
to the magnetic field, respectively. Further X{a) is the 
dephasing rate density. Note that dephasing does not 
change the energy of the spin-system. It is important to 
distinguish between this model and the two level atom 
in the previous subsection. Here the Lindblad operators 
are defined in the basis of the time-dependent eigenstates 
of the Hamiltonian whereas before the Lindblad opera- 
tors have been defined in the basis of the excited and the 
ground state of the two level atom which are independent 
of the properties of the applied laserfield and hence not 
the eigenstates of the Hamiltonian. Such dephasing oper- 
ators could be realized by random fluctuations of the field 
strength of the applied magnetic field. Since the Hamil- 
tonian changes in time, the dephasing operators have to 
be time-dependent, too. As in Section ^ we define the 
operator 



A{tMO) =\eit)) 
A{t)\gm =1.9(0) 

and find, from Eq. jnj, for p^ = A^ (t) p{t) A{t) 



(54) 



/w -i 



E 
~E 



1 



2|q gia )p (0 (q g-^a 



2p\t) 



X A (a) da 







-2iE ^A 



put) 



Wp^iit) 



(e- 

Je--l)pfi(i) 

X A(a) da. 



i)pUt) 







(55) 



Here pf^ [t) denotes the components of p^ [t) . The solu- 
tion of Eq. ^^ is 

Pn(i)=Pn(0) (56) 

pf2(i)=pf2(0)e-^(^+^)'e-/* (57) 



with / and g defined in Eqs. H47() and H48() . respec- 
tively. As the last step we need to calculate p{t) — 
A{t)p'^{t)A'^ [t). If the evolution of the Hamiltonian is 
cyclic we have A{t) =diag(e*''',e~*''') where ip is the geo- 
metric phase for |e) . The geometric phase is half of the 
solid angle enclosed by the path which |e(t)) drives on 
the Bloch sphere P| . This is equivalent to half of the solid 
angle enclosed by the path determined by the direction 
of the magnetic field. Hence we finally get for the com- 
ponents of the density operator after the Hamiltonian 
undergoes a closed loop 



Pii(r)-pn(0) 



(58) 



In the latter equation we can see a phase change due to 
the energy difference of the system, an additional phase 
change due to the dephasing and the geometric phase. 
Furthermore we see how the absolute value of the off- 
diagonal element of the density operator decreases expo- 
nentially in time because of the dephasing. 

Our task now is to remove the dynamic phase. We 
do a cra;-transformation in our system and then in the 
time interval [T, 2T] we drive the direction of the mag- 
netic field around the same loop as before but backwards: 
B{T + t) = B{T-t). The components of the density op- 
erator p'(T) after the cra;-transformation are 

p'niT) ^ P22{T) = I - pn{0) 

MT) = P2.{T) = pt.(0)e<^+^^-Me-^/ (59) 

When we drive the magnetic field backwards, then the 
parallel transported eigenstates are |e(r-fi)) — \e{T—t)) 
and|g(r + i)) ^\g{T-t)). 
Now we define the operator 



A'{T + t)\e{T))^\e{T + t)) = \e{T-t)) 
A\T + t)\g{T))^\g{T + t)) = \g{T-t)) 



(60) 



which parallel transports the eigenstates of the Hamilto- 
nian H{T + t). Again we transform the density operator 
p"^{T + t) = A"! {T + t)p'{T + t) A' {T + t) and find ioT the 
components of p"^{2T) Eq. ifSSjl 

p[i{2T) = p[,iT) = l~pn{0) 

p[i{2T) = p'i2(T)e-<^+3)^e-^/ = p,2{0re~''^e-^^f 

(61) 

where in the last step Eq. H59|l is used. Now we need 
to find p'{2T) = A'{2T)p'^{2T)A"'{2T). From Eq. ^, 
and because now we drive the loop backwards and hence 
get the same geometric phase up to a sign, it follows that 
A'{2T) = At(T) =diag(e-*'^,e+*'^) 

After another CTaj-transformation we finally get the den- 
sity operator 



p{2T) 



( Pii(O) 



Pi2(0)e4*^e-2r/^ 
l-Pii(O) , 



(62) 



Hence we see that not only the dynamic phase of the 
Hamiltonian is removed, but also the phase shift through 
dephasing. What stays is twice the difference of the geo- 
metric phases of the ground state and the excited state. 
This geometric effect appears in the off-diagonal compo- 
nents of the density operator and is damped out expo- 
nentially in time through dephasing. 



IV. GENERALIZATIONS 
A. Non-Cyclic Evolution 

To consider a non-cyclic evolution we first outline 
Pati's analysis 0|. If H{T) = H{0), Pati compared 
the phase of the parallel transported eigenstate of the 
Hamiltonian at time t = T, \n(T)) with the phase of 
the eigenstate at time t = 0, |n(0)) . If the Hamilto- 
nian does not undergo a closed loop, i.e. H{T) ^ H{0), 
then |n(T)) is not |?^(0)) up to a geometric phase. Com- 
paring the phases of states which differ not only by a 
phase is not straightforward. Pati introduced a reference 
section \n{t)) which is supported by eigenstates of the 
Hamiltonian H{t). The phase of \n{t)) is fixed by the re- 
quirement to make \n{t)) in phase with |n(0)) as defined 
by means of the work of 3, i.e. (n(0)|n(t)) = 0. Then 
\n{T)) and \n{T)) differ only by a phase and this phase 
is defined to be the generalization of the geometric phase 
to non-cylic evolution. 

We use this idea and generalize it to open systems. 
First one has to calculate the density operator p^{T) in 
the rotating axis representation as in Section^ Instead 
of transforming back to the original picture we transform 
to the picture given by the reference section introduced 
in [ij. We define the operator A(T) by 

|n(T)) ^A{T)\nm (63) 

The density operator in this new picture is 

/(T) ^A\T)A{T)p^{T)A\T)A{T) 

=diag(e*'^\ . . . , e*'^")p'^(T)diag(e-"^\ . . . , e"*"^") 

(64) 

and A'^{T)A{T) = diag(e^'^i , . . . , e^'^") is the generalized 
holonomy transformation with respect to the reference 
section \n{T)) . 

B. Non-Abelian Holonomies 

Again we consider the master equation Eq. l(Tn)l with 
an adiabatically changing Hamiltonian. For simplicity 
we restrict to a cyclic Hamiltonian H{T) — H{0). 

Until now we assumed the eigenenergies of the Hamil- 
tonian to be non-degenerate. However, if the eigenvalues 
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are degenerate we expect to get non-Abelian holonomies 
[a, 13 as the generalization of the geometric phase. In 
this case we have the eigenvalue equation 



H{t)\n^{t)) ^ Er,{t)\n„,{t)) 



(65) 



in which n — 1, ■ ■ ■ , N ; m — 1, ■ ■ ■ , Mn and M„ is the 
degree of degeneracy of the subspace of the Hamiltonian 
with energy En- The M„ are required to be constant in 
time, i.e. we do not allow any level crossings of the Hamil- 
tonian H{t). The \nm{t)) are now assumed to satisfy the 
modified parallel transport condition [8J 

(n„(i)|-|n,„,(i)) =0 Vm,m' = l,--- ,M„. (66) 
Now we can define the operator A{t) by 

A{t)\nr,M) = Wr,m (67) 

and with its help we transform the master equation in 
the rotating axis representation to remove the time de- 
pendence of the eigenspaces of the Hamiltonian analo- 
gous to Section ^ Doing this we get the new master 
equation Eq. lO for p^{t) = A^{t)p{t)A{t) as before. 
Again it can be shown that the terms p'^ {t) A'^ [t) A{t) 
and A'^ {t) A{t) p"^ {t) can be neglected if the |rim(i)) sat- 
isfy the modified parallel transport condition Eq. (|66|) in 
the adiabatic and weak damping limit. This justifies the 
condition Eq. H66|) . Since the proof for this is much the 
same as for non-degenerate Hamiltonians as done in the 
Appendix, we do not carry out the proof in this case. 
Now we have to solve (16) which represents the dynam- 
ics. Finally we have to transform back to the original 
picture: 



p{T) = A{T)p^{T)A\T). 



(68) 



We obtain the non-Abelian holonomy A{T) e U{Mi) (g) 
■■■(E) U{Mpf) similar to the way obtained the geometric 
phase for non-degenerate systems. 



V. CONCLUSIONS AND FUTURE 
DIRECTIONS 

We have introduced the rotating axis transformation, 
in which parallel transport of the eigenstates of the 
Hamiltonian plays an important role, to study the geo- 
metric phase for an adiabatically evolving multilevel sys- 
tem. This transformation was shown to be particularly 
useful in simplifying the calculation of open-system evo- 
lution, described by a master equation of the Lindblad 
form, as it allows an easy separation of dynamic and ge- 
ometric phases. These advantages were illustrated by 
applying it to optical resonance with spontaneous emis- 
sion, where we obtain known results but more easily. The 
method was then used to quickly and easily study the ef- 
fects of the geometric phase in a number of new problems. 



In one application we show explicitly how to remove the 
dynamic phase. 

Although, in our applications, we concentrated on 
Abelian holonomies for nondegenerate systems, the gen- 
eralization to non-Abelian holonomies for degenerate 
Hamiltonians and to non-cyclic evolution is straight for- 
ward. 
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APPENDIX A: NEGLECTING A''{t)A{t) IN THE 
ADIABATIC APPROXIMATION 

Here we prove that the terms containing A^{t)A{t) in 
Eq. I|12|) can be neglected in the adiabatic approximation. 
The proof will be analogous to the proof of the adiabatic 
theorem given in [27] . For simplicity we assume H{t) in 
Eq. H12|) to be diagonal which can always be achieved by 
a proper time independent transformation and hence is 
no restriction. We start by transforming Eq. H12|) in the 
interaction picture. We define 

(A^Afit) = e-^.''o ^"Wd*At(t)i(t)e^/o ^"(*)'^* 

and get by use of Eq. (|12|) the master equation in the 
interaction picture 

p" ^ p^iA^A)" - (A^Afp" + i Y: ^rs ip"] (Al) 

Q = l 

The formal solution of this equation is 



p^it) ^p"io) + j^ ip^iAur (A^Arp" 

+ ^J2^rs[p"]\^)ds. (A2) 



Within the integral, there are some contributions of the 
form of a product of a slowly varying function and a 
fast oscillating function. These contributions are known 
to become small when the frequency of the oscillating 
function increases in comparison with the time derivative 
of the slowly varying function. To see this we make use 



of the result, following [2 



Furthermore we can rewrite Eq. HA4p as 



' — > (A3) 

To make use of this we write Eq. (jA2|l in components 



4(0 - 4(0) + j^ [p^ ^^^^% - i^^^)" p" 

k 
+ 2 Z^ ^^ a»fe /^fei ^ al] ^ aik ^ akl Pi] 



pSrf^/r«J(.)d.. (A4) 



Here and later, summations are implied over all indices 
except of i and j. The components of A and T are as- 
sumed to be small in comparison with ujij = Ei — Ej 
(adiabaticity and weak damping, respectively) and hence 
we see in Eq. (|A4ll that all components of p^ are small 
in comparison with Uij. The components of (A^A) are 
slowly varying and hence the off-diagonal components, 
{A^ A)^j, i ^ j are oscillating with frequency uJij and can 
be neglected. Because of the parallel transport condition 
the diagonal elements {A'^ A)^^ are null as we can see: 

(^ ^)nn ~(^ A)nn 

= (n(0)|Ati|n(0)) = {n{t)\j^\n{t)) ^0 

The last equality is true because we assumed the \n{t)) 
to be parallel transported. So we have proved that we 
can neglect A^^A in Eq. ((T^ . 



4(0=4(0)+ / (i 



fe 

a=l 



^J- aik Pkl J- ajl 



— J- ^1., i 



^H H 



otki akl rlj 



p^'ErS;rfy)(s)ds. (A5) 



The star denotes complex conjugation. The functions 
r^r^* are oscillating with frequency w^fe ^^ji- Hence we 
can achieve a significant simplification if the differences 
of all eigenfrequences, Uik — loji are not vanishing (are 
big in comparison with H and r^) which is always the 
case if we consider a two- level system. Then the condition 
LOik—ojji = implies i = k,j = l or i=j,k = l and hence 
only corresponding parts will contribute to the integral: 



p"At) ^p"Ao) + 




■^"ij^ aikPkk^ aik 



25- V" o" t"* + 2r'f^ ■ o" r" 

^^n^ null, yto ^ -T-11 \^ ■'-'-*- r^tt /-'■■n ^ . 



-^■^J Oiii rii ^ aii 



ail rij OL'jj 



(A6) 



Here we see that only the absolute value of the off- 
diagonal elements of T^ and hence F^ contribute in 
Eq. ^. 

We can further note that if we set i = j, we find that 
the diagonal elements of the density operator are coupled 
only to diagonal elements, whereas for i ^ j, we find that 
p^j is coupled only to itself. 
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